Measurement of the 25* Hyperfine Interval in Atomic Hydrogen 
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An optical measurement of the 2S hyperfine interval in atomic hydrogen using two-photon spec- 
troscopy of the 15'-25' transition gives a value of 177 556 834.3(6.7) Hz. The uncertainty is 2.4 times 
smaller than achieved by our group in 2003 and more than 4 times smaller than for any indepen- 
dent radio- frequency measurement. The specific combination of the 25* and IS hyperfine intervals 
predicted by QED theory D21 = 8/hfs(25') — /hfs(1*S') = 48 953(3) Hz is in good agreement with 
the value of 48 923(54) Hz obtained from this experiment. 

PACS numbers: 12.20.Fv, 32.10.Fn, 32.30.Jc, 42.62.Fi 



Precision spectroscopy in simple atomic systems and 
predictions by quantum-electrodynamics theory (QED) 
supply essential data for determination of fundamental 
constants (P,!!!,!^)- The measurements also enable sensi- 
tive tests of the theory ensuring high confidence in calcu- 
lations. One type of these tests is based on the hyperfine 
splitting (HFS) measurements in hydrogen-like systems. 

In conventional atoms the sensitivity of the HFS-based 
tests is restricted by an insufficient knowledge of their nu- 
clear structure [4] . The corresponding uncertainty can be 
reduced by construction of the specific difference of the 
2S' and IS HFS frequencies 1^21 = 8/hfs(25') -/hfs(15') 
for which the nuclear size effects significantly cancel out 
[1,0, [7|]. The interest in the D21 calculations was inspired 
by experiments of P. Kusch et al ^] and significant 
progress has been made recently ^ . The present theoret- 
ical uncertainty is due to fourth-order QED corrections 
such as a{Za)'^m/M o?(Zo?) in units of /hfs(I'S') {a is 
the fine structure constant, Z is the nuclear charge, m/M 
is the electron-to-nucleus mass ratio) which are related 
to HFS tests in muonium, spectroscopy in positronium 
and the Lamb shift in H [5]. 

The IS and 2S HFS frequencies are accurately mea- 
sured in H, D, and ^He+. The lowest D21 relative un- 
certainty of 10 ppb (normalized by /hfs(15')) is reached 
in the He+ ion M while for H and D the uncer- 
tainty is about 100 ppb. The present sensitivity of these 
QED tests is mostly restricted by the experimental un- 
certainty of /hfs(25'), since /hfs(15') is measured with 
much higher accuracy. 

In 2003 we implemented an optical method for mea- 
suring the 25* HFS frequency in H by two-photon spec- 
troscopy of the 1S-2S transition pjj]. The result of 
177 556 860(16) Hz improved the previous values [1, [13] 
measured by radio-frequency spectroscopy. The calcu- 
lation of D21 from the two most recent measurements 
lH, [1^ and the precisely measured value of /hfs(I^) = 
1420 405 751.768(1) Hz [14] showed a deviation from the 
theoretical prediction (L)^f = 48.953(3) kHz, ^) at the 
2(7 level (cr is the D21 uncertainty). 

In this work we have re- measured /hfs(25') in H. The 
optical method relies on the measurement of the fre- 



quency difference between the two two-photon transi- 
tions IS{F = 0) ^ 2S{F = 0) (the singlet, /s) and 
1S{F = 1) ^ 2S{F = 1) (the triplet, ft) recorded se- 
quentially in time [12]. The splitting is obtained from 
/hfs(25') = /hfs(15') + /t - /s in zero magnetic field. 

As the most critical improvement compared to |12| w e 
use a new ultra-stable optical frequency reference [l5|. 
Also, a re-analysis of the 2S HFS frequency pressure shift 
shows that it is negligible in our apparatus. 

To measure /hfs(25') we sequentially excite the sin- 
glet and the triplet transitions by the second harmonic of 
a 486 nm dye laser [16] locked to the Ultra-Low Expan- 
sion (ULE) glass reference cavity 1 with horizontal axis 
(Figllj). The dye laser has a line width of 60 Hz (for 0.2 s) 
and a frequency drift of about 1 Hz/s while its frequency 
stability is 5 x 10~^^ for 10^ s (linear drift corrected). 

To tune the dye laser frequency between the two hy- 
perfine transitions, a double pass acousto-optic modula- 
tor (AOM) is installed between cavity 1 and the laser. 
The required frequency detuning of 310 MHz is too big 
to tune the laser rapidly without taking it out of lock. 

We take advantage of the excellent frequency stabil- 
ity of an external cavity diode laser (ECDL) at 972 nm 
locked to the thermally- and vibrationally-compensated 
ULE cavity 2 described in details in [l5[. It is stabi- 
lized at the zero expansion temperature such that the in- 
fluence of external temperature fluctuations are strongly 
suppressed. The frequency drift is nearly linear with a 
slope of about +50 mHz/s. The ECDL has a line width of 
0.5 Hz and a frequency stability of 4 x 10~^^ in 10^ s. The 
ECDL is continuously locked to cavity 2. Its frequency 
is monitored by a fiber frequency comb referenced to an 
active H-maser to ensure correct operation. 

The beatnote frequency between the dye laser and the 
second harmonic of the ECDL is mixed down with a lo- 
cal oscillator (186 MHz) to the 30 MHz range and then 
counted with a counter as shown in FigHJ All oscillators 
are referenced to the GPS-disciplined H-maser guarante- 
ing a fractional frequency uncertainty of 10"^"^. 

The measurement of the beatnote frequency shows 
that the re- locking of the dye laser causes two effects: 
First, after re- locking, the dye laser frequency shows a 
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FIG. 1: Experimental setup. The dye laser locked to the ULE 
cavity 1 in horizontal configuration serves for spectroscopy 
of the 15'-25' transition in H. To switch between the singlet 
and the triplet transitions the double-pass AOM frequency 
is changed. The 972 nm ECDL is continuously locked to the 
ULE cavity 2 in vertical configuration and serves as an optical 
frequency reference. The beat note frequencies of 32 MHz 
(triplet) and 29 MHz (singlet) are recorded by the counter 
simultaneously with the 15'-25' data. 



strong non-linear drift up to a few Hz/s due to transient 
thermal effects in cavity 1. Second, it causes abrupt ran- 
dom changes of the laser frequency of up to 100 Hz at 
486 nm due to insufficient servo loop amplification at low 
frequencies. In 2003 we had no possibility to monitor the 
dye laser frequency during the measurement. The laser 
instability resulted in excessive data scatter and could 
cause a systematic shift. In the new measurement the 
instability of the dye laser has basically no influence on 
the data quality, as the ECDL is never unlocked. 

The second harmonic of the dye laser is coupled to a 
linear enhancement cavity forming a standing wave for 
Doppler-free two-photon spectroscopy. Atomic hydrogen 
produced in a microwave discharge is cooled down to 4- 
7K by a copper nozzle and escapes along the cavity axis. 
Some of the atoms are excited to the 25* state during their 
flight. The 2S state is quenched in the detection region at 
a distance of 19.5 cm from the nozzle by an electric field of 
10 V/cm (quench field 1). The 2S population is measured 
by counting the Lyman-a photons. Compensation coils 
together with /i-metal shielding suppress magnetic fields 
in the interrogation volume to less than 30 mG. A part of 
the volume (18 cm in length) is surrounded by additional 
magnetic shielding with a suppression factor of 300 which 
also serves as beam collimator and Faraday cage. This 
shielding separates the high- vacuum region from the rest 
of the vacuum chamber. All parts of the apparatus sur- 



rounding the beam are coated with graphite to suppress 
stray electric fields. 

Compared to [HI we use a cryogenic nozzle with bigger 
diameter of 2.2 mm which allows recording H lines up 
to 1 hour without melting the frozen layer of H2 on its 
inner surface. The film improves the 25* count rate, but 
at some thickness it blocks the 243 nm laser beam and 
should be melted. An electrode is installed close to the 
nozzle to optionally quench the 2S atoms excited within 
that region (quench field 2). 

We use time-of-flight detection to study different ve- 
locity groups in the cold atomic beam. The 243 nm 
light is chopped at 160 Hz, and a multi-channel scaler 
records counts falling into 12 time bins starting at r = 
10, 210 , 2210 /is and finishing at 3 ms after the 243 nm 
light is closed by the chopper. 

The UV power in the enhancement cavity is main- 
tained nearly constant at the level of 300 mW per direc- 
tion is monitored by a calibrated photodiode installed af- 
ter the output coupler. To prevent beam-pointing effects 
we use an integrating sphere. The power fluctuations 
are taken into account by correction of the measured fre- 
quencies for the AC Stark shift which is evaluated by a 
Monte-Carlo code [I^. For r = 810 /is the AC Stark shift 
correction equals 1.3(1) Hz/mW. 

The measurement sequence is the same as in |12| : 
Groups of 2-4 singlet or triplet spectra are recorded one 
after another, the time to and the AOM frequency corre- 
sponding to each line center are defined by a fit. We use 
either a Lorentzian fit or an unsymmetrical line shape 
obtained by averaging all superimposed and amplitude- 
normalized experimental line shapes for each r ("aver- 
aged line fit"). The lines are fitted by the averaged line 
fit for the given r using 3 parameters: the amplitude 
A, the frequency offset /o and the background. The dif- 
ference ft — fs is obtained from a double linear fit of 4 
neighboring groups of values fo{to) after the AC Stark 
shift correction. A constant linear drift of the ECDL 
frequency is assumed during recording of each 4 groups. 

During 17 days of measurement in February- April 2008 
about 1200 1S-2S hydrogen spectra have been recorded 
in 28 sets (Fig. [2]). Three types of tests were performed: 
(i) variation of the particle flow coming to the cold noz- 
zle in the range (0.8-10) X 10^^ s-\ (ii) the quench field 2 
switched on/off, (iii) the direction of the compensation 
magnetic field is reversed. Further we consider the most 
important systematic effects. 

The collisional shift. In 2003 we set the upper bound 
for the 25* HFS interval frequency shift as the total shift 
of the 2S level (8MHz/mbar, [lH). Now we theoret- 
ically re-analyze the shift which appears only in third 
order of perturbation theory if collisions with H(IS') are 
considered ^18|]. For each of the excited P-states of the 
colliding partners a and h (m and n correspondingly) the 
contribution to the shift of the 2S hyperfine component 
scales inversely to the sum ~ ^mP + -^is ~ -^nP • 
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FIG. 2: 2S HFS frequency in H evaluated for r = 810 /iS 
using the averaged line fit. Each data point is the average 
over one set of measurements at unchanged conditions (par- 
ticle flow, quench field 2, compensation magnetic field). The 
histogram shows /hfs(25') data for low and moderate par- 
ticle flows (gray bars) and all data (white bars) except the 
ones taken at inverted compensation field. The solid line is a 
Gaussian fit. 
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FIG. 3: Top: Averaged experimental /hfs(25') results for dif- 
ferent r (squares). Only low and moderate flow data from 
Fig. [2] are analyzed. The gray la area shows the expected fre- 
quency shift for AT = 0.4 K. Inset: The measured amplitude 
ratio of the triplet and singlet transitions which is consistent 
with AT = 0.4 K. Bottom: mutual deviation of averaged HFS 
interval frequencies for different evaluation methods. Error 
bars are similar to the ones from the upper plot. 



differential HFS shift is on the order of 10"'^ of the 26' 
level shift due to the small ratio of the 2S HFS energy 
and the difference E^^ — E^p. The short-range inter- 
action does not contribute since for impact parameters 
smaller than 20ao (ao is the Bohr radius) the 2S state 
quenches. Integration over the discrete spectrum gives a 
result on the order of lOHz/mbar, the continuum gives 
a similar contribution. The effect from collisions with 
H2 should be of the same size because the closest dipole 
transition from the H2 ground state lies is the UV region. 
The intra-beam pressure in the nozzle is 10~^ mbar and 
rapidly decreases in the expanding beam so we can ne- 
glect the shift. A stronger effect may be induced by the 
dipole interaction with photoionized H Since there 
are only 10 protons present in the excitation zone at a 
time the effect is negligible. 

We extrapolate the data of Fig. [2] to zero flow sepa- 
rately for each of the time bins. No systematic deviation 
between the extrapolated and the averaged values is ob- 
served. For r = 810 /is the difference equals 1 Hz with 
uncertainties of 11 Hz and 6 Hz. The data taken at higher 
flows have an excessive scatter due to an instability of the 
overloaded cryogenic vacuum system, so we exclude them 
(shifting the final value by +0.5<j) from the analysis. The 
data are averaged without adding any systematics. 

Line shape / beam temperature. We analyze the aver- 
aged low- and moderate-flow data (Fig. [2]) for different 
delays r (Figl3l top) which shows an increase of the mea- 
sured /hfs(2S') frequency for shorter delays. The differ- 
ential 2nd order Doppler effect is on the mHz level and 
cannot explain the difference. The analysis of the line 



shapes indicates that the 2S{F = 0) atoms have a higher 
temperature T than the 2S{F = 1) ones (assuming a 
Maxwellian distribution). The amplitude ratio of the two 
hyperfine transitions varies from 3.1 to 3.6 depending on 
r (see the inset) which means that the fraction of slow 
singlet atoms is less than for the triplet ones. 

Fitting the singlet and triplet lines by line shapes simu- 
lated for different beam temperatures [17] we find for the 
difference AT = (Ts) - (Tt) = 0.4 K with (Ts) = 4.2 K. 
The thermalization of hydrogen atoms on the H2 film 
depends on its spin state which we still cannot explain. 

We have evaluated the expected shift of /hfs(25') 
vs. AT for different r by Monte- Carlo simulations. 
The dependencies are nearly linear with slopes of 
220(40) Hz/K (10 /is), 68(5) Hz/K (210 /is), 17(4) Hz/K 
(410 /is), 5(3) Hz/K (610 /is), and 0(2) Hz/K for higher 
delays (Fig.[3|). The uncertainties result from numerical 
errors and possible jitter of r at a level of 20 jas. As ex- 
pected, the effect vanishes at higher r since the velocity 
distribution of the delayed atoms is insensitive to the ini- 
tial distribution. The optimal compromise between the 
statistical uncertainty and the effect of different beam 
temperatures is reached at r = 810 /is. 

The influence of the fitting model and data weight- 
ing is analyzed in Fig. [3] (bottom). The Lorentzian fit 
of strongly asymmetrical spectra for r = 10 /is results in 
a shift of 7(10) Hz. On the other hand, the line shape 
for r = 810 /is is indistinguishable from a Lorentzian, so 
the effect reduces to the sub-hertz level. For r = 810 /is 
we get /hfs(25') of 177 556 835.3(6.2) Hz and add 2 Hz 
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FIG. 4: Experimental and theoretical values for the D21 
8/hfs(25') — /hfs(1*S') difference in atomic hydrogen. 



uncertainty for line shape / beam temperature shifts. 

The A C Stark shift The differential AC Stark shift of 
the 25* hyperfine components is negligible (~ 1 /iHz/mW 
p^). Though the measured difference in singlet and 
triplet excitation powers is taken into account for each 
individual spectrum, the uncertainty of correction itself 
may cause an error in /hfs(25'). Comparison of the 
data evaluated with/without correction shows a system- 
atic difference from 2 Hz to 5 Hz depending on the delay 
r. Assuming an error in the AC Stark shift evaluation 
of 30% (including the error in the power measurement), 
the contribution to the final uncertainty budget is 1.3 Hz. 

The DC Stark shift. We have no possibility to mea- 
sure stray electric fields in our apparatus and evaluate its 
contribution of —1(1) Hz as in [l2|. The influence of the 
quench fields 1,2 is tested with the help of simulations 
[l7|, an effect on the sub- hertz level is expected. 

Magnetic fields. The sensitivity of /hfs(25') to an ex- 
ternal magnetic field B equals +9 600 5^ Hz/G^. In 5 
sets of measurements the compensation field direction 
was inverted which increased the field in the less shielded 
zone from 30 mC to 300 mC (Fig.[2j). The corresponding 
/hfs(25') differs from the value measured with proper 
orientation by —3(12) Hz at r = 810 /is. We estimate the 
uncertainty resulting from magnetic fields as 0.5 Hz. 

Summarizing the uncertainties (Table[Ij) we get the fi- 
nal result /hfs(26') = 177 556 834.3(6.7) Hz. The corre- 
sponding value D21 = 48 923(54) Hz is in good agreement 
with the theoretical prediction (Fig.[4|) which is the first 
result after consistent with theory within 1 a. This 
measurement is part of the long-term project on preci- 
sion spectroscopy of the 15'-2S' transition in H [l6| and 
has an impact on its centroid frequency used for the Ry- 



dberg constant derivation [2[. 
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TABLE I: Uncertainty budget for the new 2S HFS frequency 
measurement in atomic hydrogen. 





Frequency [Hz] 


Uncertainty [Hz] 


Averaged interval frequency 


177 556 835.3 


6.2 


Line shape/temperature 





2 


DC Stark shift 


-1 


1 


AC Stark shift 





1.3 


Magnetic fields 





0.5 


Final result 


177 556 834.3 


6.7 
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